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ABSTRACT 

Cold fronts ~ sharp discontinuities recently discovered by Chandra in many clusters of 
galaxies - are believed to be due to a hot gas flow over a colder gravitationally bound 
gas cloud. We analyze the stability of the fronts with respect to Kelvin-Helmholtz 
instability and show that an intrinsic width of the interface of the order of a few per 
cent of the curvature radius strongly limits the growth of perturbation. For the best 
studied case of a front in the Cluster Abell 3667 we conclude that current observational 
data on the width and extent of the front can be explained even in the absence of 
dynamically important magnetic fields. 
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1 INTRODUCTION 



Cold fronts were discovered as sharp features in the X-ray 
' surface brightness distribution in Chandra observations of 
; the clusters A2142 and A3667 (Markevitch et al., 2000, 
• Vikhlinin, Markevitch & Murray, 2001a), see also Marke- 
■ vitch et al. (2002). Similar features have now been found 
I in several other clusters (e.g. Sun et al., 2002, Kempner, 
Sarazin & Ricker, 2002). Unlike shocks, these features have 
lower gas temperature on the X-ray brighter side of the dis- 
continuity. For that reason they are called "cold fronts". It 
is believed that some cold fronts are formed when a sub- 
cluster merges with another cluster and the ram pressure 
of gas flowing outside the subcluster gives the contact dis- 
continuity the characteristic curved shape. Indeed, features 
resembling cold fronts are found in numerical simulations of 
cluster formation (Bialck, Evrard & Mohr, 2002, Nagai & 
Kravtsov, 2003). Ablation of the gaseous cloud by the hot 
gas causes characteristic differential motion of the gas inside 
the subcluster, which transports the low entropy gas from 
the subcluster core towards the contact discontinuity, thus 
enhancing the jump in temperature and surface brightness 
across the discontinuity (Heinz et al. 2003). 

Here we address the question of the front stability. As 
was pointed out by Vikhlinin et al. (2001a, b), Vikhlinin and 
Markevitch (2002) the observed fronts are narrow and could 
be unstable to the Kelvin-Helmholtz (KH) instability. For 
A3667 however the front appears to be narrow (less than ~5 
kpc in width) up to ~30°from the stagnation point. Mag- 
netic fields can act as a stabilizing agent thus allowing indi- 
rect estimates of the field strength near the front (Vikhlinin, 



Markevitch & Murray,2001b). On the other hand numerical 
simulations without magnetic field (and in particular rela- 
tively high resolution simulations by Heinz et al. 2003) do 
not show instability of the front within 20-30°. While this 
discrepancy could be due to numerical effects, we revisit the 
question of the front stability below and show that the char- 
acteristic convex geometry and finite (small) intrinsic width 
of the interface may help to stabilize the discontinuity with 
respect to KH instability. 

The structure of the paper is as follows. In Section 2 
we derive a simple expression for the growth of the KH in- 
stability along the curved interface. In Section 3 we argue 
that diffusion processes are likely to set an approximately 
constant width for the interface. In Section 4 we discuss as- 
trophysical applications. The last section summarizes our 
findings. 



2 KH INSTABILITY ALONG THE 
SPHERICAL INTERFACE 

For an infinitely thin plane parallel interface separating 
semi-infinite layers of two incompressible fiuids with den- 
sities pi and p2 the dispersion relation for Kelvin-Helmholtz 
instability is (Landau & Lifshitz, 1987): 



ui = kv 



1 + ' 



(1) 



where A: = 27r/A is a wave number, A is the wavelength of the 
perturbation, v is the velocity of fiuid with density pi, while 
fiuid with density p2 is at rest and M = The frequency cu 
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Figure 1. Geometry of the problem. Hot gas is flowing around 
an approximately spherical eloud of colder gas. Stream lines illus- 
trate the hot gas pattern of motion. The gas inside the cloud is 
assumed to be still. While the interface is unstable against KHI, 
observations shows that the front remains thin up to angles of 
~30°. Perturbation with a wave number ki at angle Bi will have 
a wave number k2 when arriving at angle 62 due to the increase 
of the velocity along the interface. 



is measured in the frame of the fluid at rest (i.e. with density 
p2). For cold fronts in clusters we adopt the geometry of the 
flow following Vikhlinin et al. 2001a as shown in Fig.l. The 
front has a curvature radius R and the velocity v along the 
interface is assumed to follow the law v = vo sin 9, where 
6 is the angle along the interface. For the potential flow of 
incompressible fluid vo = 1.5 Voo, where Voo is the velocity 
of the cold gas through the hot gas. For a particular case of 
the cold front in A3667, Vikhlinin et el. (2001b) argue that 
due to compressibility vo is more close to Voo ■ For simplicity 
we assume that the densities of both fluids remain constant 
along the interface. 

Consider the perturbation of the interface separating 
cold and hot gas. From eq.(l) it follows that perturbations 
are advected along the interface with the increasing velocity 



Va = VaosinO, where Va 



Vo 



1+tl 



This impUes that the 



wave number k of the perturbation decreases as: 



sm fi 



(2) 



Va sm o 

where fci = k{9i) is the wave number at some initial posi- 
tion 01. This result can be formally obtained from the low- 
est order WKB (Wentzel-Kramer-Brillouin) approximation 
(see appendix), applicable if kRsinO ^ 1. Apart from the 
exponentially growing factor, the next order WKB approx- 
imation allows one to estimate slower changes of the per- 
turbation amplitude d caused by the changes of the velocity 
along the interface and stretching of fluid elements in the az- 
imuthal direction (see appendix). For estimates we assume 



that <5 oc sin 6'-'", and a ^ 2. I.e. 5 = di where 5i 

is the initial amplitude. 

Thus one can evaluate the growth of the perturbation 
propagating along the interface as: 



_ /sin(9iV 



exp 



7(t)dt } , 



(3) 



where j{t) is the increment of KH instability and one has 
to evaluate it using k from eq.(2). An analogous expression 
was used by Inogamov & Chekhlov (1991) for the growth 
of perturbations for Ray leigh- Taylor instability. From eq. 
(1) and (2) it is clear that the increment 7 oc few does not 
increase as the perturbation propagates, since the increase 
of the velocity amplitude is compensated by the decrease of 
the wave number. Thus 'y{t) = feiUfosin^i = const, where 
Vio = voj^. Replacing dt with dt — j^dd = ^dO one can 
write an explicit expression for J '){t)dt: 



j{t)dt = 



R 



tiVio sm Oi- 



Vao sm f 



Rki sin6ii— ln[tan(0/2)] . 



(4) 



Thus the final expression for the growth factor for a pertur- 
bation starting at 61 with the wave number fei measured at 
position 02 is: 



/ amei \ 
Vsin^y' 



• r 1 

exp < Jifci sin 61 — — 



Growth factor = 
ln[tan(e/2)]|^J 



(5) 



One can compare this expression with eq. (10) from 
Vikhlinin & Markevitch (2002) - i) an additional factor ap- 
pears in front of the exponential and more importantly ii) 
the argument of the exponential is changed due to the vari- 
ation of the velocity along the interface, thus reducing the 
growth factor for a fixed fei and 61. We note here that 
from the point of view of the degree of nonlinearity of a 
perturbation one can use the ratio of the amplitude to the 
wavelength: S/\ — S-^. Therefore the expression for the 
degree of nonlinearity of a perturbation will contain an ad- 
ditional factor proportional to ^ = in front of the 
exponential. 

Alternatively, since fcisin^i — k^ sin 62, we may regard 
the growth factor as a function of k2 specified at the angle 
02- The formal condition for the growth factor (5) to di- 
verge for small 0i is Rk2 sm02/^/]I > 2. Thus in spite of the 
suppression of the growth factor for a fixed initial wave 
number, the growth factor can be infinitely large for a fixed 
final wave number. There are however several factors which 
can limit the growth of instability. We argue below that a 
small, but finite, thickness of the interface is the most im- 
portant factor in the conditions relevant for cold fronts in 
clusters. 



2.1 Finite thickness of the interface 

Any real interface can not be infinitely thin. If the thickness 

of the interface is h, then only modes with k < kmax ^ 
1/h are unstable. Thus for a fixed 02 and fc2 < kmax the 
minimum angle at which KH instability starts is such that 



fc2 



sin 02- In Section 3 we argue that difiiusive 
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Wavelength A/R 

Figure 2. Growth factor as a function of the wavelength A2 for 
several values of angle 62 and intrinsic with h. X and h arc ex- 
pressed in units of the front curvature radius. The factor in 
eq.(7) is set to 1.43. 

processes across the interface may set a finite thickness of 
the interface, which is almost independent of the position 
(angle). One can approximately account for the finite width 
of the interface by using Omin as a lower bound in eq.(5). But 
in fact the largest contribution to the growth factor comes 
from small angles and this part of the integral has to be 
evaluated more accurately than simply introducing a cutoff 
in 9. One can use for instance the dependence of increment 
on the wave number and thickness in the form given by 
Rayleigh (see "The Theory of Sound", e.g. 1945 edition) for 
two fluids with similar densities: 

l{k, h) = ^ ^e-'^^h _ (fe/j _ 1)2^ (6) 

where k < 1.2785//i is the condition for instability. Thus the 
growth factor is: 

GF=f^)'exp(:^i^ r ^e-^^^^-ikh-l)^d0\. 

\ Sm 62/ [h VaO Jg^ J 

where 9i = 6rmn is such that k2^^ = 1.2785//i. In this 
case the growth factor has to be calculated numerically. In 
Fig.2 we show the total growth factor calculated for several 
values of 02 and h. For 02 ~ 30°the instability grows very 
strongly if the thickness of the interface is less than --^1.5-2% 
of the curvature radius. We note that all the above expres- 
sions for the growth factor do not depend explicitly on the 
absolute value of vo , indicating that the appearance of the 
front may be similar for clouds moving with different veloc- 
ities, if the relation between the width and the curvature of 
the interface is the same. 

An interesting question is if the nonlinear evolution of 
KH instability itself can be responsible for establishing the 



width of the interface (along the lines of reasoning given by 
Nulsen, 1982) of order of few percent of R, smoothing an 
interface and suppressing growth of longer modes. We leave 
this question for subsequent studies and in the next section 
consider the width of the interface set by diffusion. 



3 WIDTH OF THE FRONT DUE TO 
TRANSPORT PROCESSES 

For flow past a flat plate (with high Reynolds number) the 
width of the boundary layer is known to vary as the square 
root of the distance x from the front edge of the plate (Bla/- 

sius law): 

h - yJxvjU, (8) 

where v is the viscosity. For a rounded body, the flow velocity 
increases linearly with distance from the stagnation point, 
U ~ vqx/R^ making the width of the boundary layer con- 
stant (Landau & Lifshitz, 1987). Thus h ~ sJvRJvq, where 
_R is the curvature radius. Substituting further v ^ ViXi, 
where Vi and A, are the rms velocity and mean free path 
for ions and assuming that Vi ~ vo one gets an order of 
magnitude estimate of the expected width h ^ \/\iR. In 
the absence of magnetic field the viscosity of the hot gas is 
larger than for the cold gas and most of the velocity drop 
may take place on the hotter side in the form of a broad 
layer. We note here that observationally it is almost impos- 
sible to "see" such a layer. 

Analogously one can consider the case of a scalar diffu- 
sion (e.g. thermal conduction) with a constant diffusion coef- 
ficient D. Assuming that conditions inside the cloud (where 
the fluid is not moving) are frozen and seeking a stationary 
solution of a diffusion equation in outer gas it is easy to show 

that h - .[^. 

Both expressions, when applied to the case of A3667 and 
compared with the upper limit on the front width (Vikhlinin 
et al., 2001a), indicate significant suppression of transport 
processes across the front. More accurate calculations and a 
detailed comparison with observations are beyond the scope 
of this paper and will be reported elsewhere. In what follows 
we assume that the interface has a constant width indepen- 
dent of and the expression (7) for the growth factor is 
applicable. 



We now apply the above results to the particular case of 
the cold front iu A,"5667. We set the following parameters 
based on the results of Vikhlinin et al., 2001a, b, Vikhlinin 
& Markevitch 2002: R = 410 kpc, 6I2 = 30°, ^ = 1.43. 
They give the width of the front (upper limit) in terms of 
the standard deviation of a gaussian smoothing, as ct < 5 
kpc. In our notation, h — \/27ra ~ 2.5(t. In Fig. 3 we plot 
the growth factor calculated for several values of a and angle 
62- This is largely the same figure as Fig.2 but now in natu- 
ral units. For comparison we show (upper dotted curve) the 
growth factor calculated for an infinitely thin interface, as- 
suming fixed wave number k and O2 =30° . One can see that 
accounting for varying wave number and the flnite thickness 



(7) 

4 APPLICATION TO A3667 
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4.2 Role of gravity and compressibility 
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Figure 3. Growth factor for A3667 parameters. Observational 
limit on the width of the interface is cr < 5kpc. Upper dotted 
curve corresponds to the estimate of the growth rate for fixed k. 



of the interface (within the observational limits, i.e. a < 5 
kpc) drastically reduces the growth factor. For these reduced 
values it is not obvious that small initial perturbations can 
grow strongly nonlinearly and disturb the sharp appearance 
of the front around 62 ~30° . For larger angles or for smaller 
widths of the interface the growth factors are large. 



4.1 Role of the cloud internal motions 

While we assumed above that the gas inside the cloud is 
still, internal motions are likely to appear during the for- 
mation of the front and especially during later periods of 
the cloud evolution when continuous stripping at the sides 
of the cloud is compensated by the cold gas flow along the 
interface to replenish the losses. E.g. in simulations by Heinz 
et al. (2003) the velocity of the gas inside the cloud ~0.2 of 
the velocity of the surrounding gas at the moment when the 
morphology of the front and the temperature distribution 
resemble those observed in A3667^. The internal motions 
have a two-fold effect: i) they reduce the shear rate and ii) 
they increase the group speed of the perturbation in the 
interface, both of which reduce the total growth of a per- 
turbation arriving at the angle 62- Accounting for internal 
motions can easily reduce the argument of the exponential 
in equation (7) by 30-50% thus strongly reducing the final 
growth factor compared to Fig.3. 



^ We note here that in the absence of gravity the velocities of the 
colder and hotter fluids would be related as Vcoid = 'Vhot^/V' ^ 
follows from the Bernoulli equation. Due to gravity the velocity 
of the motions in cold gas should be lower. 



In the above discussion we have completely neglected the 
role of gravity, although we assumed that it plays an impor- 
tant role in preserving the integrity of the cold cloud when it 
passes through the hot gas. For an infinitely thin interface, 
gravity modifies the dispersion relation (1) to: 



kv- 



1 



kl + ^l (1 + m)' 



(9) 



where g is gravitational acceleration. For a fixed k2 and 62 
and given that k = sin 62 / sin 9i and v — vo sin the role 
of gravity is more important for small angles 6, effectively 
introducing a lower limit on Oi . The growth factor calculated 
with the above dispersion relation and using the values of 
g and vo from Vikhlinin & Markevitch (2002) is very large, 
implying that gravity does not strongly suppress KH instar 
bility, in agreement with their conclusion. 

For an interface of finite thickness, gravity makes the 
region near the stagnation point stable for perturbations of 
any k. One can estimate the Richardson number of the flow 
assuming that density and velocity both vary over a layer of 
thickness h: 



Ri ' 



-a dp 
p dz 

w 



Apgh 

n 7)2 



(10) 



The sufficient condition for stability is Ri > 0.25 everywhere 
in the flow. Again plugging in the numbers from Vikhlinin & 
Markevitch (2002) we flnd that for cr = 1, 3, 5 kpc the front 
is stable up to angles of ~6,ll,14°respectively. Simple esti- 
mates show that account for this stable region would further 
reduce the growth factor compared to Fig.3, although not 
dramatically. This is expected given that gravity alone (for 
the thin interface case - sec above) does not strongly sup- 
press KH instability. Accurate calculations of the growth fac- 
tor taking account of gravity and the finite thickness of the 
interface require knowledge of the density and velocity pro- 
files of the interface and are beyond the scope of this paper. 
We note only that stretching of the flow in the stable region 
near the stagnation point insures that initial perturbations 
in the unstable regions further downstream are small. 

Since we are considering motion near the stagnation 
point the velocities are small compared to the gas sound 
velocity (even for colder gas) and compressibility docs not 
play a significant role. Further downstream the compressibil- 
ity starts to be more important and it reduces the growth 
factor of the perturbation (see dispersion relation in Miles, 
1958, Gerwin, 1968 or Nulsen 1982). E.g. for the density ra- 
tio n = 0.5 the quantity is equal to = 1.41 for 
small shear velocity and it drops to ~ 1.13 when the shear 
velocity approaches the sound speed of the colder gas^ . 



5 CONCLUSIONS 

We have shown that accounting for wave number changes 
along the curved interface and finite (but small) width of 
the interface significantly reduces the KH instability growth 



^ An argument of the exponential in eq. (5) is proportional to this 
factor 
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factor. We also argue that if diffusion sets tfie intrinsic width 
of the interface then this width does not vary much along 
the interface (at least for small angles). 

For a set of parameters relevant for the front in A3667 
the growth factor is not large enough to guarantee nonlinear 
growth at angles ~30°, provided that the intrinsic width of 
the front is of the order of a few percents of the curvature 
radius (i.e. not much less than the existing observational 
limit). Therefore dynamically important magnetic fields may 
not be necessary to stabilize this particular cold front. The 
growth factor (and therefore possible limits on the magnetic 
fields) depends critically on the width of the interface 
and much more weakly on the angular extent over which 
the front remains sharp. A factor of 2-3 improvement in the 
observational constraints on the interface width is therefore 
extremely important. 

For a given width of the interface the value of the growth 
factor does not depend strongly on the absolute value of the 
shear velocity. This probably explains the ubiquity of fronts 
both in observations and hydrodynamical simulations. 
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APPENDIX A: WKB APPROXIMATION FOR 
PLANE ACCELERATING FLOWS 

Consider a plane boundary between two incompressible, in- 
viscid, irrotational fiuids with densities pi and p2 moving 
in the x direction with velocities ui{x,y) and U2{x,y). The 
fluid '1' is accelerating along the interface: = ui/x = 
const. To match pressures at the boundary we have to as- 
sume that the velocity of fluid '2' is U2 — ui ^/pi/p2 = Ui^/j2 
and is also accelerating along the interface. We now consider 
small perturbations of the boundary £,{x,t) and velocity po- 
tentials of two fluids <pi^2{x,y,t). The governing equations 
are: 



A01 



: 



901 



dy 



+ Ux,l,2^ 4 — ^ 

ot ox ay 



(Al) 
(A2) 



(A3) 



902 d4>2 

dt ' dx 

where Ux.i:2 and ttj/.i,2 arc the velocity components in x 
and y directions for two flows. The second equation is the 
kinematic boundary condition and the third equation is the 
Bernoulli equation. These equations are evaluated at the 
boundary ^. 

We are seeking a solution of these equations 
in the form £, = 6(a;)e^(*'' (^»/^ and 0i,2 = 
Ol,2(a;,^/)e*(^*W+^*'"^)+■^''•l•=^''^ which decays for y = ±oo, 
where a, 6, and all S axe slowly varying functions and e 
is a small parameter. Using these expressions and keeping 
only the terms of order 1/e wc recover a dispersion relation, 
which we write in a "spatial" form, denoting ^§f-/£ = —OJ, 

plUx,l + P2Ux,2 ± i{Ux,l — Ux,2)^/PlP2 



PiuIa + P2ulo 



(A4) 



Now collecting higher order terms and using the relations 
g^'^ = ±i-|ip from the Laplace equation and = 
piMx.i aoj^ from the Bernoulli equation we get a transport 

P2Uj. 2 OX ^ o X- 

equation for spatial variations of the amplitude b: 

I 2 2 \ db , dUy 1 dUy 2 . , „ / J _N 

(PlMx.l + '°2Mx,2)^ - KPlUx,!—^^ + P2Ux,2—^)b = 0. (A5) 

Thus 



dx 



(A6) 



Incomprcssibility implies that = — = —Ux/x. The 
last equality is due to assumed form of acceleration. There- 
fore 



b oc X 



(A7) 



i.e. the amplitude decreases with the distance as 1/x. If ad- 
ditional uniform stretching in the 3rd dimension is present 
then one has to set = — — . In the case of inter- 
est for us (potential flow past a sphere) one can simply set 



dy 

b (X x~ 



-2ux/x. Substituting this into eq. (A6) we get 



(A8) 



